Classical Rotons in Cold Atomic Traps 
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We predict the emergence of a roton minimum in the dispersion relation of elementary excitations 
in cold atomic gases in the presence of diffusive light. In large magneto-topical traps, multiple- 
scattering of light is responsible for the collective behavior of the system, which is associated to 
an effective Coulomb-like interaction between the atoms. In optically thick clouds, the re-scattered 
light undergoes diffusive propagation, which is responsible for a stochastic short-range force acting 
on the atoms. We show that the dynamical competition between these two forces results on a new 
polariton mode, which exhibits a roton minimum. Making use of Feynman's formula for the static 
structure factor, we show that the roton minimum is related to the appearance of long-range order 
in the system. 
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Since the first ideas advanced by Landau [2 [5], the 
concept of the "roton minimum" in the dispersion of the 
collective modes of a certain physical system has played a 
central role in the description of superfluidity. After the 
success of the theory in the context of superfluid phases 
of 4 He, rotons have received considerable attention since 
then and have been identified in many different quantum 
interacting systems. Recently, Cormack et al. [3] suggest 
that rotons may appear in already moderately interact- 
ing ultracold Bose gases and Kalman et al. [I] have nu- 
merically observed their emergence in two-dimensional 
dipolar bosonic gases. 

In fact, the emergence of a roton minimum in the ex- 
citation spectrum strongly depends upon the shape of 
the interacting potential or, equivalently, on how parti- 
cles are correlated. The correlational origin of the roton 
minimum has been firstly suggested by Feynman, where 
the static structure factor S(k) is expressed in terms of 
the dispersion relation ui(k) as 
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This has an enormous implication on the interpretation 
of the physical properties of the system in terms of the 
dispersion relation: the presence of a roton minimum is 
the signature of strong correlations system. In the limit 
of large mode softening, i.e., for rotons with zero fre- 
quency, the system can develop mechanical instabilities, 
which can lead to interesting physics phenomena. In Ref. 
[5] , Henkel et al have suggested that the presence of a "ro- 
ton zero" is at the origin of crystallization in ultracold 
Rydberg gases. 

In this letter, we describe the classical origin of a ro- 
ton minimum in the excitation spectrum of cold atomic 
clouds confined in magneto-optical traps (MOTs). Due 
to competition between long-range interactions between 
the atoms and the stochastic forces associated to the dif- 



fusion of light, atoms in MOTs experience a complex ef- 
fective interaction, which we find to be associated with a 
polariton dispersion relation. This polariton mode is re- 
sult of the dynamical coupling of the density waves with 
the fluctuations of the light intensity inside the trap. 

A route for the most intriguing complex behavior in 
large magneto-optical traps relies exactly in the multiple 
scattering of light [BJ [7] . Due to the consecutive scatter- 
ing and re-absorption of photons, the atoms experience 
a mediated long-range interaction potential similar to 
Coulomb system (~ 1/r) 019] and the system can there- 
fore be regarded as a one-component trapped plasma. In 
a series of previous works, we have put in evidence the 
important consequences of such plasma description of a 
cold atoms gas [TOl E] , whereas the formal analogy and 
the application of plasma physics techniques reveal to be 
important in the description of driven mechanical insta- 
bilities |12H16j or even more exciting instability phenom- 
ena, like phonon lasing [17J . Moreover, in such optically 
thick traps, it is known that the light does not propagate 
ballistically, rather exhibiting a diffusive behavior |18j . 
In this situation, the energy transport velocity Ue, i.e. 
the velocity that accounts for the propagation of energy 
by the scattered wave, is smaller than c [19, 20]. Labeyrie 
et al. |2jQ have experimentally observed that Ve can, in- 
deed, be several orders of magnitude smaller than c in the 
case of resonant light propagating in traps, already with a 
moderate optical thickness, thus putting in evidence the 
phenomenon of slow light. More recently, the diffusive 
behavior of light has been identified as a source of dy- 
namical instabilities leading to the formation of photon 
bubbles in magneto-optical traps [2"2"] . 

In what follows, we consider that the dynamics of cold 
atoms in MOTs is described by the Vlasov equation 



dt 



V+i)]F r V v /(r,M)=G, (2) 



where f(r,v,t) is the normalized distribution function 
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1= dr dvf(r,v,t). 



(3) 
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The total force J2i Fi = Ft + F c accounts for both the 
trapping and cooling forces. There are evidences [23T 
|2"5] that the density profile is approximately constant for 
large traps (typically with N ~ 10 9 — 10 10 atoms), which 
allows us to consider the system to be homogeneous and 
thus to neglect the effects of the trap. The collective force 
can be described by a Poisson equation [91 [TU] 



V • F c (r,i) = Q ct! J dvf(r,v,t), 



(4) 



The pre-factor in Eq. Q represents an effective charge 
Qdt = &l{&r — &l)I/c of the atoms induced by light, 
where or. and ox represent the scattering and absorp- 
tion cross sections [6 9J, and / is the light intensity. For 
most of the experimental conditions, the scattering cross 
section is larger than the absorption cross section, i.e. 
cr > enforcing the effective charge to be a positive 
quantity. We have showed that the positiveness of Q e{{ is 
an essential condition for the existence of stable oscilla- 
tions in the system (see e.g. Ref. [TU]). 

We now assume that the diffusive behavior of light 
inside the trap can be macroscopically described by the 
diffusion equation 



SI-- 
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(10) 



where 
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is the photon inhomogeneity parameter. Combining the 
latter result with Eq. ([7|, we finally obtain the kinetic 
dispersion relation 



1 = 



m - v k 2 



i 



dfo 



— Q/k dv x 



dv, (12) 



where we have considered perturbations parallel to the 
wave- vector k = fce z , for definiteness. Here, we have 
defined two typical frequencies of the system. The first 
one is associated with the oscillations of the atoms due to 
the long-range force, corresponding to an effective plasma 
frequency [10] 



dt 



V -VVI = 0. 



(5) 



The diffusion coefficient is determined by T> = £ 2 /t , 
where the photon mean free pass is I — with 
n = no J fdv standing for the atomic density. According 
to experimental results |21j . the diffusion time r can be 
considered as independent from the atom density, so the 
diffusion coefficient explicitly reads 



UJ,, 



(13) 



The second important quantity is the rate at which the 
photons scatter inside the trap, or simply the diffusion 
frequency 



f3n 2V 2 / 



■D Q . 



(14) 



V(r,t) 



T 2 ^-2^2 



J /(r,v,t)dv 



(6) 



We now linearize the Eqs. ([2|, Q and ^ by allowing 
fluctuations around the equilibrium values / = /o + Sf, 
I = Iq + SI and V = V Q + SV, such that 



d_ 
at 



v • V I Sf + -SF C ■ V v / = 



(7) 



We notice that this frequency depends on the scale at 
which the diffusive processes occur (micro-, meso- or 
macroscopic), as it depends upon the spatial scale L 
at which the light intensity varies. We will discuss the 
macroscopic case below. 

The integral in Eq. (12 1 can be evaluated using the 



Landau prescription, according to which the full informa- 
tion about the initial conditions is cast if the integration 
path is set to pass below the pole f2 = v z k. We split the 
integral into two parts 
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(8) 



(9) 



where Q oft ,o = o-l{o~r — <jl)Iq/c. Assuming peri- 
odic perturbations on both space and time, such that 
(Sf, 51, 5V) cx exp(ik • r - iClt), Eqs. (|6) and ^ yield 



i dfo dv = pr r i dfo dv 

v z — f2/ k dv z J v z — Vt/k dv z 



+J7T 



dfo 

dv z 



-n/k 



(15) 



where Pr stands for the Cauchy principal value. Assum- 
ing a phase speed w ph = uj/k much greater than the width 
of the distribution, such that fo and its derivatives get 
small as v z gets large, we may expand the denominator 
in ( 15 ) which, together with the relation 
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FIG. 1: (Color online) Illustration of the real (top panel) and 
imaginary (bottom panel) parts of the polariton dispersion re- 
lation in the macroscopic regime for 2?o = 2.0\j J Lo p . We can 
observe the emergence of a roton minimum for u>d = 1.9cj p 
(full black line) and aid = 1.99cj p . The rotons softens the 
frequency to uj(k Iot ) = at the critical value u) d c ^ = 2.0uj p . 
Roton instability is illustrated here for cj = 2.2ui v (red line). 
The short-wavelength oscillations corresponding to A: Ad > 1 
are kinematically damped. The usual plasma dispersion rela- 
tion of ref. [TU] is presented here for comparison (dashed gray 
line). 
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dfo 



v z — fl/k dv z 
simply yields 



fo 



( Vz - n/kf 



(16) 



1 df 
v z — fl/k dv z 



k 2 v 2 \ 

fo ( 1 + J dv x dv y dv z . 

(17) 

Assuming the atomic equilibrium to be described by a 
Maxwell distribution 



(18) 



with v th = yJksT ' /m standing for the thermal speed, we 
may finally write 



n 2 



u 2 s k 2 



iil - Vnk 2 



+iir 
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k 2 dv? 



(19) 



-n/k 



where we have defined the sound atomic speed u s — 
V3v th . Separating the frequency into its real and imag- 
inary parts, n = uj + «7, with 7 <C w, we may finally 
write 



= (u 2 p + u 2 s k 2 ) 1 



u d V k 2 



w2 + W 



(20) 



and 



7 



, ,2 

UJ d u p 



w| + V 2 k± A; 3 A 



,-3/(2fe 2 A^) 



(21) 



where Xd = Ms/w p is the effective Debye length. This 
dispersion relation describes a quasi-particle excitation 
resulting from the atom-photon coupling, or a polariton. 
This is the main result of this paper and we now explicitly 
show that is contains a roton minimum. 

The diffusive description of light is the result of a 
macroscopic approximation which is known to hold if ab- 
sorption takes place at scales much larger than the mean 
free path £ |18j . This is true provided the following hier- 
archy for the relevant length scales 



A < i < a < L, 



(22) 



where A is the light wavelength and a is the size of the sys- 
tem. According to typical experimental conditions |21j . 
the mean-free path is found to value t ~ 300 fim and 
the diffusion coefficient T> ~ 0.66 m 2 s _1 . Based on our 
previous estimates |10j . the effective plasma frequency 
and Debye length respectively value u> p ~ 27r x 100 Hz 
and Xd ~ 100 fim. Therefore, provided the identifica- 
tion L = 2V 2 I /I in (14| and using the inequalities 
in Eq. (p2l with a 



1 mm 



the diffusive approx- 
imation is valid if the diffusion and plasma frequencies 
are of the same order, i.e. uid ~ w p . This is attained if 
Id = Dq j uj p ~ 1 cm is of the same order of the intensity 
variation length L, which is a reasonable condition for 
typical experimental scenarios. Nevertheless, by varying 
the number of atoms in the trap it is possible to control 
the optical thickness of the system and, therefore, tune 
the value of the diffusion coefficient T> , which makes our 
present estimates even more flexible. In extremis, it may 
be also possible to attain different diffusion regimes, but 
this situation is out of the scope of the present paper. 
In Fig. 0, it is shown that a roton minimum emerges 
in the dispersion relation ( 20 ) in the diffusive regime. 
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As the values of ojd increase (i.e., for stronger diffusion), 
the frequency decreases around the roton wavenumber 
k Iot ~ Xd/^j,- This feature is often referred to as mode 

softening. At the critical value oj^ — 2w p , the mode soft- 
ens towards zero, which is a clear manifestation of a roton 

fc) 

instability mechanism. For Wd > oj \ , the system enters a 
crystallization phase. This mechanism has been recently 
discussed in the literature as it can lead to the forma- 
tion of supersolids [5 . An important remark is related 
to the Landau damping at short wavelengths. Modes in 
the region k\u ^ 1 undergo a kinematic damping. For- 
tunately, rotons are possible to be excited at longer wave- 
lengths (k rot \[) < 1), thus avoiding the Landau damping 
mechanism. Moreover, the onset of diffusion tends to de- 
crease the damping rate (see Fig. ([I])). This nourishes 
hope for rotons to be experimentally observable. 



A remarkable feature of the polariton spectrum in ( 20 ) 



is that it exhibits a roton minimum in a three dimensional 
system even in the absence of strong interactions. It is 
clear from the application of Landau's criterion that the 
present spectrum does not correspond to that of a su- 
perfluid, as the mode is gapped at the origin. This is 
a consequence of the long-range nature of the interac- 
tion between nearby atoms, and the quasi-particle mass 
is proportional to uj p . This therefore corresponds to an 
important example where rotons are not intrinsically con- 
nected to Goldstone modes. 

Another important property of the classical rotons de- 
scribed above is that they carry useful information about 
the long-range correlation of the system. By using the 
extension of Feynman's formula ([I]) to finite temperature 
systems [57], the static structure factor S(k) is given by 



S(k) 



hk 2 



2mu>(k) 



coth 



/ Hu)(k)\ u 
\2k B T 



uj{kf 



(23) 



where we have used the non-degeneracy condition fc^T 3> 
huj. In this limit, the structure factor corresponds to that 
obtained based on a hydrodynamic treatment [28] . In 
Fig. ([2]), we illustrate the behavior of S(k) for the same 
parameters of Fig. (JlJ . 

The static two-point correlation function g(r) = 
(n(r)n(0))/(n(r))(n(0)) ~ (n(r)n(r' ')) / 'n 2 , can then be 
easily calculated provided the relation g(r) = 1 + 
J 7-1 [S(k) — 1] 29J, which after the integrating out the 
angular variables simply reads 



7(r) 



1 



k sin(fcr) 



[S(k) - 1] . (24) 



As it can be observed in Fig. the appearance of a 
minimum in the excitation spectrum ( 20 1 is associated 



with the occurrence of long-range correlation in the sys- 
tem. By inspection, one founds that the correlation func- 
tion oscillates with the period T = 2ir/k [ot . This feature 
can be qualitatively understood in the context of Percus- 
Yevick theory [301 EI] , where the correlation function is 




FIG. 2: Static structure factor S(k) depicted for Do 



1.9uj p (black full 



2.0\jjujp. uid = (gray dashed line), ujd 1 
line) and uid = 1.99o;p (black dashed line). A pick emerges in 
the static structure factor around k Iot — Xu/i\ = O.SA^ 1 . 



approximated by g(r) FV ~ 1 + cor _1 cos(fco + So)e~ K ° r , 
where cq is a constant and zq = KQ + iko is the pole of the 
function S(k) — 1. We remark, however, that the PY the- 
ory was originally developed for hard-sphere potentials, 
and therefore does not describe systems with long-range 
interactions. For that reason, we have not used it to 
compute g(r). 
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3: Two-point correlation function g(r) depicted for Do = 
(gray dashed line), uid = 1.9cj p (black full 
line) and u)d = 1.99cj p (black dashed line). The correlation 
function oscillates with the period a of 27r/fc rot ~ . 
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In conclusion, we have derived the dispersion relation 
for the atom-photon polariton in large magneto-optical 
traps in the presence of diffusive light. We have explic- 
itly computed the excitation spectrum for the particu- 
lar case of a thermal atomic distribution, revealing the 
emergence of a roton minimum for a set of parameters 
compatible with current experimental conditions. We 
have also shown that the increase of the light diffusiv- 
ity lead to mode softening around the roton minimum, 
which may eventually drive the system into a roton in- 
stability situation. Using the relation between the static 
structure factor and the dispersion relation, we have ex- 
plicitly demonstrated that the roton minimum is related 
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to the emergence of long-range correlations in the system. number SFRH/BD/37452/2007. 
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